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I. INTRODUCTION 

The study of the properties of strongly interacting 
quark-gluon plasma (QGP) 1 is one of the most inves- 
tigated research topics in recent years. The interest is fu- 
eled on the one hand by direct experimental questions re- 
lated to the properties of QGP produced at RHIC and on 
the other hand by the possibility of studying analytically the 
non-perturbative properties of plasmas in various, mostly su- 
persymmetric, gauge theories using the AdS/CFT correspon- 
dence [2 . Although there does not exist so far a direct coun- 
terpart of the QCD plasma, the study of exact properties of 
similar systems not based on ad-hoc phenomenological mod- 
els may increase our understanding of the properties of the 
QCD plasma. This hope has been substantiated, e.g., by the 
discovery of a universal strong coupling shear viscosity to en- 
tropy ratio, which is valid for a wide range of different gauge 
theories [3,4. 

A lot of work has been initially undertaken for a static 
plasma system with a fixed constant temperature. After the 
earliest investigations of transport coefficients various other 
observables have been discussed related to the properties of 
fundamental flavors in the thermal medium such as drag force 
calculations [5, 6 , jet-quenching 7, 8, 9 , meson spectrum 
TO, 11 , meson melting 12 , thermodynamics of fundamental 
flavors 13 . 

On the other hand, since the experimentally produced 
plasma is always a non-static expanding system, it was tempt- 
ing to extend the AdS/CFT investigations to such a time- 
dependent dynamical setting. Qualitative duals for ther- 
malization and cooling have been suggested in 14, 15 , a 
quantitative framework for studying boost-invariant expan- 
sion of a plasma system in — 4 SYM has been pro- 
posed in 16 . Subsequent work within this framework in- 
clude 17 , 18 , 191 1^ 1^ ESI m [251 [26], which concentrated 
on the details of the dynamics of the expanding plasma in 
pure A" = 4 SYM theory. 

The main motivation for this work was to bring together the 
two lines of investigation and to study an expanding plasma 
in the A/" = 4 SYM theory with fundamental flavors. The 
flavor system is represented by an embedding of a D7 brane 
in the dual geometry of the expanding plasma, which is - in 
an essential way - time-dependent. Such time-dependent em- 
beddings have not been considered so far in the literature on 
flavor systems. In this paper we would like to concentrate 
on the features of the system that are the direct analogues 
of the corresponding properties studied in the static case. In 
particular we find the (time-dependent) embedding, we find 



the time-dependence of the chiral condensate including lead- 
ing viscosity effects and we describe the behavior of mesonic 
modes in the scalar and vector channels. Finally we perform 
a holographic renormalization of the D7 action density. We 
will show that the chiral condensate, meson spectra and ac- 
tion density are compatible with the adiabatic approximation, 
i.e. their leading contribution agrees with the result obtained 
from the static AdS/Blackhole by naively assuming Bjorken 
scaling. 

As a word of caution let us emphasize that the gauge theory 
that we are considering does not exhibit spontaneous chiral 
symmetry breaking so from this point of view differs substan- 
tially from QCD. However, since this is the first investigation 
of a flavor system in an expanding plasma, we chose to deal 
with the simplest theoretical setting, which is the only sys- 
tem for which the dual expanding geometry is known so far. 
Ultimately one would like to extend these investigations to 
theories that exhibit chiral symmetry breaking. 

The plan of this paper is as follows. We will first re- 
view some basic facts about viscous hydrodynamics from field 
theory perspective and its implementation in AdS/CFT. In 
the next section we will perturbatively determine the time- 
dependent embedding of a D7 brane in that geometry and 
determine the chiral quark condensate. Thereafter, we deter- 
mine meson spectra from fiuctuations about the embedding. 
Moreover, we give the regularized D7 action and show that it 
is again compatible with expectations from adiabatic consid- 
erations. 



A. Boost invariant kinematics 

An interesting kinematical regime of the expanding plasma 
is the so-called central rapidity region. There, as was sug- 
gested by Bjorken 27 , one assumes that the system is in- 
variant under longitudinal boosts. This assumption is in fact 
commonly used in realistic hydrodynamic simulations of QGP 
28 . If in addition we assume no dependence on transverse 
coordinates (a limit of infinitely large nuclei) the dynamics 
simplifies enormously. 

In order to study boost-invariant plasma configurations it 
is convenient to pass from Minkowski coordinates (x", x^,x±) 
to proper-time/spacetime rapidity ones {T,y,x±) through 

x^^rcosht/ a;^=rsinhj/ (1) 

The object of this work is to describe the spacetime depen- 
dence of the energy-momentum tensor of a boost-invariant 
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plasma in A/" = 4 SYM theory at strong coupling. The sym- 
metries of the problem reduce the number of independent 
components of T^^ to three. Energy-momentum conservation 
9^1 = and tracelessness = allows to express all 
components in terms of just a single function - the energy 
density e(r) in the local rest frame. Explicitly we have 



T^T = e(r) 



T 



(2) 
(3) 

(4) 



Gauge theory dynamics should now pick out a definite func- 
tion £{t). 



B. Viscous hydrodynamics 

The object of a hydrodynamic model is to determine the 
spacetime dependence of the energy-momentum tensor for an 
expanding (plasma) system. The simplest dynamical assump- 
tion is that of a perfect fluid. This amounts to assuming that 
the energy momentum has the form 

Tfj,^ = (e + p)u^u^ + pr)^„ (5) 

where is the local 4- velocity of the fluid {v? = —1), £ 
is the energy density and p is the pressure. In the case of 
N = A SYM theory that we consider here T/^* = and hence 
e = 3p. The equation of motion that one obtains from energy 
conservation in the boost-invariant setup is 



drE — 



4 e 



(6) 



C. AdS/CFT description of an expanding 
boost-invariant plasma 

In 16, 17, 18, 19 the program of perturbatively determin- 
ing the dual geometry to the viscous hydrodynamic model 
discussed in the previous section was carried out. The non- 
compact part of the metric takes the form 



L2 



1 dz^ 
= \~9^.dX^ dX" + ^ 



(11) 



with z the holographic direction. The energy-momentum ten- 
sor is related to the fourth order term of expanded in z 
32 . 



1 

— lim — r( 
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(12) 



Of course it is too difficult to perform this construction 
for arbitrary functions e(r). What has been performed in 
practice is an expansion of e(r) for (large) proper-times r 
16, 17, 18, 19 . 

Each subsequent term in the expansion can then be deter- 
mined by requiring regularity of the square of the Riemann 
tensor order-by-order. In it was shown that in order to 
study the large proper-time limit of the metric one is led to 
introduce a scaling variable 



(13) 



and take the limit r — > oo with ii fixed. For the discussion 
of subleading terms in the metric an expansion around this 
limit was performed, 



whose solution is the celebrated Bjorken result 



A{z,t) = 00(1;) -I- a\{v) — 



i 2 



+ . 



(14) 



£0 
4 

r 3 



(7) 



and similarly for the other coefficients. The leading and first 
subleading coefficients are given by 



Once one wants to include dissipative effects coming from 
shear viscosity, the description becomes more complex. In 
a first approximation one adds to the perfect fiuid tensor a 
dissipative contribution r7(V^M^-|-VyU^), where is the shear 
viscosity of the fluid. The resulting equations of motion get 
modifled to 



9t£ = 



4 £ Arj 
'St'^ 3^2 



(8) 



Note that in the above equation the shear viscosity is gener- 
ically temperature dependent {rj oc T"^ in the Af = i case) 
and hence r dependent. In order to have a closed system of 
equations we have to incorporate this dependence through 



?7 = 7)0 ■ £ 



(9) 



with ?7o being some numerical coefficient. The resulting en- 
ergy density satisfying (8) gets modified from (7) to 



ao(«)=ln(i^^^, ai(^;) = 2r,o^5+^, (15) 

bo{v) = ln(l + «V3), bi{v) = -2no^ + 2r;oln |^, 

co(«) = ln(l + //3), c-i{v) = -2f,oj^ - Vo In |^. 

rjo is an undetermined integration constant (which has the 
physical interpretation as the coefficient of shear viscosity) . It 
can be fixed from non-singularity of the metric; the resulting 
value 

1 

Vo = (16) 
2234 

is in accord with the known viscosity coefficient of A/" = 4 
SYM, rj/s — l/47r, calculated in the static case [3 . 

From the leading order terms in (15), using the similarity 
with the static black hole metric, we may read off the temper- 
ature from the position of the horizon. As discussed in 18], 
to this order, the result 



£0 
4 
r 3 



1 - 



2vo 



+ . 



(10) 



T(r) = 



4£o 



1 



Vo 



(17) 
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is compatible with the Stefan-Boltzmann law. This can be 
explicitly checked by extracting the energy density from the 



metric expanded in the scaling limit up to 0[r 
(12). We obtain 



sir) 



27r2 ■ A 



2r,Q 

i 2 

En T3 



through 



(18) 



which is consistent with viscous hydrodynamic evolution |(To| 



II. TIME-DEPENDENT D7 EMBEDDING 
A. D7 embeddings 

Conventional AdS/CFT describes A/" = 4 SYM theory, 
which has only fields in the adjoint representation, since all 
fields sit in the same supermultiplet containing the gauge field. 
A standard method 30 for the introduction of quenched mat- 
ter into the correspondence is by adding probe D7 branes. 
Strings stretching between the probe and the original D3 stack 
generating the AdSs x background give rise to an A/" = 2 
hypermultiplet in the fundamental representation. 

Let us consider geometries Ms x with A/5 asymptotically 
AdSs and line element 



J 2 , 2 

dsio = ^"S4 + ^ 



V + p dni + {dx^f + [dx^Y (19) 



with = E(^')' = + + {Xy = 

Placing a D7 probe parallel to the first eight coordinates, 
its position and shape are described by the coordinates X^ 
and A". Manifestly, the D7 brane breaks the 5'I7(4) ~ 50(6) 
symmetry of the internal manifold to 50(4) x 5089(2) ~ 
SU{2)l X SC/(2)b X U{1)r. The 50(2) symmetry in the 
8,9-plane can be used to rotate the embedding of the D7 to 



A* = 0, 



X^ = <S>{p,r) 



(20) 



such that one scalar field is sufficient for a complete descrip- 
tion of the embedding. 

The embedding $(p, r) is determined from the D7 action, 

5d7 = M7 J d''^e-'^^detP[g]ab + Fat 

+ [ d\P[{dpCi)]e°'^~<A^d0A^, (21) 



where denotes the pull-back to the world-volume. 

While AdSs x permits constant embeddings "1? = m, 
more general geometries require a non-trivial profile $(p). To 
our knowledge the present article is the first-time that also 
time-dependence for the embedding has been considered in 
the context of fiavored AdS/CFT. 

Close to the boundary the embedding behaves as 



c 

7^2 



(22) 




L-TiT 



FIG. 1: Chiral condensate as function of the quark mass for 
the static AdS/Schwarzschild background. 



where m and c are related to the bare quark mass ruq and 
chiral condensate (O) by 13 ^ 



{O) 



2-Ka'' 
NfN, 



(23) 
(24) 



where A = g^M^c = ^ixg^N^ = L*/{2li). For the AdS/ 
Schwarzschild geometry such a VEV forms, but vanishes for 
TUq such that there is no spontaneous chiral symmetry 
breaking associated to its appearance 11 . 

The common procedure for finding the embedding is to 
derive the equations of motion from (21) and impose the re- 
quirement that the embedding have an interpretation as a 
holographic renormalization group fiow. In particular this 
means that $(p) should be one- valued as a function of the 
holographic energy scale r, which requires regularity as p — > 0. 
(This condition is not sufficient however as will be discussed 
below.) Because of the parity symmetry r t-^ — r of the under- 
lying geometry, which manifests itself as p 1— > — p invariance of 
the equation of motion, embedding solutions are either parity 
odd, known as 'black hole' solution since they hit the horizon, 
or even, i.e., dp<^ — 0, known as 'Minkowski-type' solutions. 
It is easy to show that the condition for Minkowski-type so- 
lutions corresponds to the absence of a conical defect on the 
brane 29 . 



B. AdS/Schwarzschild: Adiabatic Approximation 

Figure 1 shows the chiral condensate as a function of 
the quark mass for the AdS/Schwarzschild geometry. Since 
all quantities enter the equations of motion in units of the 
black hole mass mth (we express everything in terms of 
ttT = 8Gjv/(37r)mi,h), the leading order time-dependence of 
the chiral condensate can be determined in an adiabatic ap- 
proximation using the fact that the perfect fiuid geometry for 



^ To be more precise these are the mass of the Af = 2 hyper- 
multiplet and the vacuum expectation value of the operator 
O = ^■^p + q^^q + mqq^q. 
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constant r appears like an AdS/Schwarzschild solution with 
fixed temperature (17). 

Our numerical results provide c/T^ ~ (m/T)". a is about 
—6 at m ~ 3(7rT) and increases to —5.2 when going to m ~ 
8(7rT). In other words 



-(3-c)/3 



-8.2/3 



(25) 



However, as has been observed in 29 , the numerics deter- 
mining Figure 1 is not particularly accurate for large quark 
mass because the chiral condensate becomes small. Therefore 
it is not possible to increase the quark mass till the exponent 
saturates. 

For large quark mass, the solutions are very far from the 
black hole and become approximately constant embeddings 
as in the supersymmetric scenario 30 . This suggests the fol- 
lowing perturbative analysis 13 . For small e we seek regular 
solutions 



$(p) = m + e/(p) 



of the DBI action in AdS/Schwarzschild geometry [11| . 



(26) 



_ 3 / , ejL^nTf 



v/l + <l>'(p) (27) 



We obtain 



{L'^Tf (3m" -f 3p^m^ + p^) 
96m5 (m? + p^f 



1 f L^'-kT 



96 V 



(28) 
(29) 



Viscous fluid 



Restricted to the scalar $, the action reads for our geometry 







B 



$2 



(p2 +$2)2' 



A:= (1 



I exp 



B := f— exp 



(1-^)2 



9 - w« 



I r2 



r3 ^p2 $(p,t)2 ' 

1 NcNt 



(31) 
(32) 

(33) 
(34) 



where Td7 = 27r/(gs(27r^s) ) is the D7 tension. D,3 = 2tt-' 
and Vx = J dy d^x± are the volume of the unit three-sphere 
and spatial part of the boundary, respectively. The latter is 
of course infinite. 

Since the viscous fluid geometry discussed in the previ- 
ous section behaves similar to AdS/Schwarzschild with time- 
dependent temperature, we do not expect spontaneous sym- 
metry breaking either, though going to small quark masses 
leaves the domain of validity of the geometry and it is thus 
hard to make a definite statement. We will therefore only 
consider 'Minkowski-type' embeddings that avoid the horizon 
at the center of the geometry. 

For a given quark mass, in general regularity is only possi- 
ble for a discrete set of values for the chiral condensate. In the 
regime under consideration, since no phase transition occurs, 
we expect c = c(m) to be a one-valued function. 

The equation of motion arising from (31) is a non-linear 
partial differential equation. We will solve it perturbatively 
by a late-time expansion 



$(p,r) = m + ^f^{p)T- 



(35) 



We have written (29) in a way that emphasizes the known 
fact that the temperature can be effectively removed from the 
equations by a suitable redefinition of the quark mass and 
condensate. In the present context, such a rescaling is not 
desirable. With (17) we conclude that the adiabatic estimate 
for the perfect fluid geometry is 



(30) 



We use a fraction of 1/3 in the exponent because all exponents 
showing up in the background geometry (11) are integer mul- 
tiples of one third. The ansatz reduces the equations of mo- 
tion to the following (infinite) system of ordinary differential 
equations 



p-'d.ip'fUp)) 



8m£§ 



9(m2 +p2)5 



-4?7oeo 



1/4 



■Up) (36) 

ff i = 8 
ff i = 11 

else; provided i < 14 



We will check in the following that this is indeed the leading 
contribution for the expanding plasma geometry. 



The boundary behavior of solutions to (36) is 



Mp) 



> mi + 

p^oc p^ 



(37) 



2 For conformity with our conventions, we have replaced their wg 
by <I> and the parameter b by L^ttT, with T the temperature of 
the black hole. 



which becomes an exact solution when the inhomogeneous 
term vanishes, li = 0. The first term, rrii, contributes 
miT~^^^ to the bare quark mass. Since we do not accept a 
time-dependence of the bare parameters on physical grounds, 
we require = 0. Thus in conjunction with regularity the 
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FIG. 2: Embeddings with m — \e~^^^ at different times. For 
late times (bold), the supersymmetric embedding is approx- 
imated. The diagonal line in the first plot corresponds to 
"I>(r) = r, where the embeddings have to end. 



value of Ci is completely fixed. In particular = implies 
Ci = or /i = 0. 

To the considered order the solution is 



with 



^(p, t) = m + c- 



54m^ 



{m? + p2)3 



(38) 



(39) 



In figure 3 we show the effect of the subleading term on the 
shape of c(r). Note that we can trust this result only as long 
as this term is small compared to the leading order, 



A -1/4 -2/3 „ 1 



(40) 



We observe that the viscosity term enters the chiral con- 
densate exactly the same way as it enters the energy density 
squared. We may therefore express the chiral condensate as^ 



{o) = 



1 Nf>? 
2f67r4 m rrif 



(41) 



^ We would like remind the reader that e(r) ~ T'^N^, such that 
(O) ~ NcNfT^. 



FIG. 3: Time-dependence of the chiral condensate for adi- 
abatic approximation (monotonic curve) and with viscosity 
correction. The circle at r = (4/3)^''* ~ 1.38 is outside the 
regime of validity given by equation (40). 



This solution is however only valid for late times as each 
term in e(r) apparently diverges for r ^ 0. While at this 
stage it is not clear up to which time (if at all) our perturba- 
tive expansion converges, we may still ask for which range of 
T it is self-consistent. For this consideration, three time-scales 
are of potential importance. Firstly, the time when the so- 
lution touches the horizon, i.e., $(0,r) = 3"^/''L^eJ^^r"^/^. 
Secondly, the moment when the viscosity term in the expan- 
sion dominates the leading term in such a way that the em- 
bedding "recoils" and stops being a one- valued function of the 
holographic direction r. This happens when d?^{r,T) / dr'^ — > 
oo. Thirdly, the time r = (4/3)''/** ^ 1.38 when the chi- 
ral condensate changes its sign and by (41) would lead to an 
imaginary energy density. Before this time scale is reached, 
the regime of validity (40) is left. 

For numerical computations and plots, we will use e^^^^ 
(which is a length) as a unit to express dimensionful quan- 
tities. Figure 2 shows the m — 1/2£q^''^ solution at various 
times before break-down. (Again we would like to stress that 
the solutions may be invalid even before - by "break-down" 
we denote their having become invalid for sure.) Figure j4] 
compares the three effects. For small quark mass solutions 
are invalidated by touching the horizon, and by imaginary 
energy density for large quark mass. 



III. MESON SPECTRA 

In the D3/D7 framework, meson spectra are determined 
from regular, normalizable solutions to the equations obtained 
from linearizing the full equations of motion of the D7 brane 
about the embedding solution that describes the position and 
shape of the brane 10 . 

In the following section we distinguish between four dimen- 
sional meson modes, which carry a "4d" label and eight di- 
mensional fluctuations, which always start with a 6 followed 
by a (Greek or Latin) capital letter, e.g. 5$ or SA^ . Our 
ansatze are products of spherical harmonics y on the inter- 
nal manifold, wave forms parallel to the boundary and radial 
parts, which describe the dependence on the holographic co- 
ordinates and are denoted by 5 followed by a small letter. 
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FIG. 4: Break-down time of the solution as a function of the 
quark mass. The horizontal line indicates the time where the 
energy density becomes imaginary, the dashed flat curve is 
where embeddings become two-valued as a function of the 
holographic energy scale r and the steep curve shows where 
solutions touch the horizon, i.e., can no longer be considered 
'Minkowski- type'. The gray area marks invalid solutions. The 



circle is at m ~ 0.682ei 



-3/8 



which is the smallest mass that 



does not lead to solutions that eventually touch the horizon. 



A. Boost invariance for the AdS geometry 

Before turning to the actual holographic computation of 
meson spectra for the viscous fluid, it is insightful to investi- 
gate how boost invariance changes wave forms of scalar and 
vector mesons in the conventional setting, where there is no 
time-dependence of the meson mass. 

In four dimensional Minkowski space the massive Klein- 
Gordon equation assumes the form 



-drrdT + r ^di + 

T " 



(ci Jo(c^r) + C2Yo(a.r))e±"=^^^ 
fci = fc| -f kl 



(42) 
(43) 



with Jo and Yo Bessel functions of first kind. Since a linear 
combination of the Bessel functions will appear frequently in 
our expressions we introduce the short hand 



Cl 



ip^judr) + C2Yp{JujdT). 



(44) 



We will not explicitly denote the time-dependence of J^o aris- 
ing from the integral over r. At this stage, the integral has 
been chosen for later convenience and gives uj r for constant 
frequencies. The eigenfrequencies ui := ^y + k\ are to be 
determined in our holographic setup. We will thus assume 
fc_L = from the start to obtain the mass spectrum. 

The 4d meson field is given as the boundary value of (lin- 
ear, normalizable) fluctuations (5$, 5'^ about the embedding 
solution <I>(p) = m. 



X* = 0-F5*(p,r), = $(p) +<5<I>(p,r). 



(45) 



For the presentation of our ansatz, we will concentrate on the 
scalar mode 5$; the pseudoscalar mode 5^ will be treated 
analogously. With the following holographic ansatz 



the boundary value corresponding to quantum number I is 
defined by 



(47) 



Equation \4&\ is a natural modification of the ansatz given in 
[lOj to separate the D7 equation of motion in anti-de Sitter 
space: 



(p^ + m?)^ T 



1 1 3 1 

-drTdr + -^dpP dp + ^Ag3 



5'^{p,t) = 0. 

(48) 



The radial equation obtained after separation reads 



^dpp^dp 



(49) 



which is the well-known result of [lOj . For simplicity, we will 
only consider the lowest Kaluza-Klein mode on the internal 
S^ such that £ = 0,y° = 1. 

The requirements of regularity in the interior (p 0) and 
vanishing at the boundary (p — > oo), fix the modes completely. 
One obtains a discrete set of modes the lightest of which is 
given by 



^0 



m? + p^ 



(50) 



For a four-dimensional massive vector meson we have 

VaF"'' = hfAla- (51) 

We assume that the solutions are still plane waves in the 
2:2, 2;3-plane, i.e., Al^ = ^''(t) expiifcj_a;_L. This yields the fol- 
lowing component equations 



-dlAt" - UrAt" + ik2(dr + i) 



+k2k-iAf 



{Ar + ki)AT and (2 ■ 
id^{k2At'' + kaAl-") 



Ar = 



Equation (|52t) can be treated separately. Its solution is 



A 



0. 



(52) 

(53) 
>3) 
(54) 

(55) 



The others may be solved without loss of generality by 
turning the coordinate system such that fca = 0. Then it 
follows immediately that A2% = C2,3.?X) [1^2,3] expifcxa;x. With 
this modified ansatz and plugging (54) into (53) we obtain 



kl 
,,,2 



At" + k2hs ( 1 



2 ^ 

£^3 \ Aid 



3 =0 

and (2^3), (56) 



which can only be satisfied for UJ2 = (^23- Moreover, 

since it is a homogeneous system, ciJ23(M, fe, fe) can be de- 
termined from degeneracy of the coefficient matrix. We shall 
not reproduce the final expression, but just note that u>23 = M 
when ^2 = fca = 0, which could also have been obtained di- 
rectly from (53). We thus end up with the two solutions 



5$(p,r) =5</>(p)j^oH3' (5"' 



(46) 



^2,3 = 6,3-?^o[<^23l expik±x±. 



(57) 
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Therefore, we adapt the holographic ansatze for meson 
modes found in 110 as follows 



Type 
I 

II, 

Il2,3 



■ 5aiiy{p)TTi[uj]e 



a = 5, 6, 7; 



5Ac 

5Ay 



III SAo 



and (2 ^ 3) 



= 5 a I u{p) 



(58) 



with the respective other components set to zero. We will only 
consider modes of type II, which are the only modes dual to 
vector mesons and therefore most interesting. 

In the ansatze, the dependence in the 0,1,2,3 directions 
has been modified as compared to what can be found in 10 . 
The reason these changes are straight-forward is the following: 
The calculation of jlOJ only uses two important properties of 
the ansatze regarding derivatives in those directions. 



/^AdSAi = KrSAi, 
0, 



/G [0, 



,71 



glddaSAb 



(59) 
(60) 



where g^d ~ diag(— l,r^, 1, 1) in our case, whereas in 10 it 
was a Minkowski metric. For our ansatze, the gauge condition 
(60) is either trivially obeyed or follows from |54|. Moreover 
it can be used to turn (59) into (51). 



B. Viscous fluid geometry 



redefine the split (62). So in order to uniquely specify such 
an ansatz we have to supplement the usual regularity condi- 
tion at p = and Dirichlet boundary condition at p = oo by 
another condition which makes the split (62) unique. In this 
paper we will impose a condition on the profile of the mode 
5(f> induced on the boundary 



$4d(r) = lim p^5(j){p,T) 



(65) 



Namely we will set 



^idT 



(66) 



which provides a definition of our frequency uj4,d{r). For con- 
stant aJ4d(T) this reduces of course to the pure AdS^ result 
(|61|. 



Our motivation for the above form ( 66 ) is that it arises as a 
WKB approximation to a Klein-Gordon equation with time 
dependent mass spectra: 



□ $4 



-drTdr+T 9„ + 9i 

r " 



<1>4 



Mid{r)<l>A 



(67) 



Before coming to the actual holographic computation, we 
would like to discuss the general framework of late-time per- 
turbative expansions that we use. 

Since the viscous fiuid geometry and our D7 embeddings are 
time-dependent, we do not expect, and do not see, a separa- 
tion into a purely r dependent and p dependent factor. This 
makes the problem very difficult to tackle analytically. We 
are helped by the property that at late proper times the ge- 
ometry becomes pure AdS^ with the corresponding D7 brane 
embedding. In this limit the simplest solution looks like (50) 



To 



1 

m? -|- 



(61) 



where J?-o[(x;] is defined in equation (44). For smaller proper- 
times it is natural to treat the frequency appearing in (l61|() 
as depending on r. However as the equations do not allow 
for a separation of variables we have r dependence also in the 
remaining part: 

f[u{T)\f{p,r) (62) 

where we allow for a general T which should reduce to JTq [u)] 
for constant u. We have moreover the expansions 



lo{t) = lo + 



1 



.f{p,r) = r>{p) + ^f'\p) + ^ 
r 3 



(63) 
(64) 



Note that the above form is not uiuque. Redefining the co- 
efficients of the expansions in an appropriate way, we may 



We may separate variables by assuming a plane wave in the 
2, 3 plane and obtain uj'ld{T) = M|j(t) + fc^. (Though we will 
assume k±_ = 0, henceforth.) The remaining equation 



-d^rdr^idir) = u)ld{T)^id{r) 



(68) 



can only be solved approximately, e.g., by the WKB approx- 
imation, which gives two linearly independent solutions. 



UJAdT 



(69) 



The square root prefactor ensures that Abel's theorem is ful- 
filled, such that the Wronskian for our ansatz is 



W 



const 
t 



(70) 



as it should be for the exact solution. 

We have now all ingredients in place to actually calculate 
the meson spectrum for the time-dependent viscous fluid ge- 
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ometry. We expand the D7 action (21),* given by 
Cdbi 



+ 



-{dp Ay) 



- (I) + (quartic) 



-,1/2 



(I) 



(11) 
(III) 



+ {drX^f(.d,xy 

- 2(drX''){dpX:'){dpX^)[drX'') 

+ e-'(n) + ^(in) 



(71) 
(72) 



-2(a,x»)(apX»)(ap^2)(a,A2) 

+(9.4 J' + {d^X^fidrAyf + (9,X')'(9pA,,)' 

-2(a.x')(apX')(ap4j,)(9,A„) (73) 

to quadratic order in fluctuations 



X' 



$ + (5<I>, 

+ 5yl2, 



X^ = + <5*, 
4„ = + M„. 



(74) 



The resulting equation of motion is evaluated by performing 
a perturbative expansion in t~^^^ , 



<5<I. = CiJo(/a;(«(r)dr)^50,(p)- 
+ C2Yo(/c.W(r)dr)^50,(p)- 

00 



-j/3 



(75) 



(76) 



and analogously for the other fluctuations.^ We use the 
known asymptotic expansion of the Bessel functions and ob- 
tain schematicaUy the following equation 

[polynomial in r^^''^] cos(J tj''^' dr) 

+ [polynomial in r^^^^] sin(/ c^<"*' dr) = 0. (77) 

At any given order, the requirement that the coefficients of the 
polynomials vanish, provides a differential equation for &<i>i 
and (5<^i depending on oji (and lower order solutions). We have 
to go to order 6 before the viscosity 770 enters the equations. 
To this order, the equations for (5(/!>i and can be separated 
by choosing suitable linear combinations and yield &(i>i = 5cf>i, 
which is what is required for the WKB ansatz (69) to be 
applicable. We impose the boundary conditions 



0, 



p-O 



finite. 



(78) 



* We do not write out those quartic terms that can only produce 

terms quartic in fluctuations. 
^ The equation of motion of Ay requires a slightly modified ansatz 

given in the appendix. 
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FIG. 5: Late-time spectra for the viscous fluid geometry. The 
supersymmetric spectrum is shown as dashed horizontal lines. 
Scalar mesons are shown as continuous lines, vector modes are 
dashed. 



The first two of these conditions pick regular normalizable so- 
lutions; the last ensures that meson solutions on the boundary 
(65) satisfy the constraint (70). Consequently, the conditions 
fix two integration constants and the frequency Ui at each 
order in the perturbative expansion. The only remaining free 
constants are the overall factors ci and C2 of our ansatz (75). 

Each of the coefficient functions has to satisfy a differential 
equation that is best expressed with the substitutions 



2 



2 / 2 

y = -p /m , 

m^{{2n + 3f 



(79) 



1). 



The lowest order equation then reads 



[y(l - y)9y + (c - (a + b - l)y)dy - ab] My) = 



-n — 1, 



c = 2. 



(80) 



This is exactly the hypergeometric equation already encoun- 
tered in 10 . The boundary conditions fix n to be a non- 
negative integer, thus yielding a discrete meson spectrum 
M = iL!o{n) and the solutions in terms of (degenerate) hy- 
pergeometric functions 2F1 are 



fti =2Fi(-n- l,-n,2;y), (81) 
h2 = (l-y)^+'"2Fi(n + 2,n + 3,2n + 4;l-y). (82) 

Only hi is regular, such that 500 = hi. 

Higher orders in perturbation theory produce inhomoge- 
neous terms in the analogues of (jSOj). Since it is a linear 
ordinary equation, the solution can still be obtained in closed 
form by standard methods. However, the resulting integrals 
are hard to solve in general. Since both solutions hi,2 are 
rational functions of y and Iny, it is however easy to do so 
for definite n. For the lowest five mesons n = 0, . . . , 4 we give 
the solutions in the appendix. 
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The mass of the lowest scalar meson mode is 



, ,{<P) _ 47r 



^4^4/3^3 



1 - 



2vo 



(83) 



To the considered order, pseudoscalar modes Sip have ex- 
actly the same equations of motion and the spectrum is de- 
generate. Moreover we note that the spectrum agrees with 
the adiabatic approximation even including the viscosity cor- 
rections. The reason for this might be that the bulk metric 
coefflcients that enter the calculation for the scalar mesons 
can be expressed completely in terms of the energy density, 
whereas the components for the y, 2, 3 directions cannot. 

The vector mesons deviate slightly from the scalar modes. 
For comparison we plot the mass ratio of scalar and vector 
modes in Figure 7 The mass of the lowest vector mesons is 
given by 



2407r<lT*/3m3 



2407r<lT''/3m3 



18»70 
7r3e4 



(84) 



(85) 



This agrees with the adiabatic approximation excluding the 
viscosity term. Since the metric components that enter the 
holographic computation, Qyy and 322, agree only up to the 
viscosity terms, this deviation does not come as a surprise. 

We plot the five lowest meson modes in Figure 5 The 
leading order term gives the exact supersymmetric spectrum 
that is approached for t ^ 00. 



C. Comparison to the adiabatic approximation 

In this subsection, we will review some properties of low- 
temperature meson spectra for the static AdS black hole. 
Plugging in the time-dependence of the temperature into 
the static meson spectra, yields an estimate for the time- 
dependent spectrum, which we will refer to as adiabatic ap- 
proximation. We will assume that the temperature depen- 
dence given in terms of Poincare time t can be obtained from 
(17) by substituting r for t: 



T. 



AdS/BH 



4eo 
3 



1 

7rt3 



(86) 



Note however that we do not expect the resulting adiabatic 
meson spectra to accurately give the viscosity corrections. 
The reason for this is that even though the horizon posi- 
tion can be expressed completely in terms of the energy den- 
sity, such that the Stefan-Boltzmann law holds, the bulk 
metric and energy momentum tensor nevertheless contain 
additional viscosity terms that cannot be captured by the 
AdS/Schwarzschild geometry even when the geometry near 
the horizon and near the boundary is matched. 

In Figure 6 we plot the numerical solution (dots) of the 
static case. We calculate the asymptotic solution in a low 
temperature expansion (solid curves), which agrees with the 




FIG. 6: The plot shows the (pseudo) scalar meson spec- 
trum. For small temperatures, the supersymmetric spectrum 
(dashed) is approached. The solid lines are asymptotic T — > 
solutions, which are in good agreement with the numerical so- 
lutions (dots) for small temperatures. 



numerical calculation. Plugging the time dependence of the 
temperature (86) into this analytical approximation, we ob- 
tain the time-dependent meson spectrum in adiabatic approx- 
imation. The mass of the lowest scalar and vector modes are 
given by 



ad. 

4>,i' 



ad. 



47r 

71 

4-K 
4-K 

71 

4-K 

71 



320 m| 



£0 



3A^ 



SOtt^ ml 
320 m| 



2407r'* ml 



2to 
t3e4 



(87) 



The mass of all of these modes decreases for increasing tem- 
perature.^ Moreover, the adiabatic scalar modes completely 
agree with our result (83), whereas the vector modes only 
agree in the leading contribution; the viscosity terms differ. 
We consider the agreement of the scalars accidental in the 
sense that it is a consequence of a certain property of the ex- 
panding plasma geometry: All metric coefHcients entering the 
scalar equation of motion can be expressed purely in terms of 
the temperature, while additional viscosity terms only show 
up in those metric coefficients that end up in the equations 
for the vector modes. 

An important assumption in our calculation has been that 
the Klein-Gordon equation is obeyed by scalar particles. This 
assumption can actually be proved employing the holographic 
equation of motion resulting from the linearization procedure. 
There is a parity symmetry p ~p in the equation of motion. 
Since only 'Minkowski-type' solutions are considered, this will 



defined by equation ||66[l 



A similar decrease of meson masses for increasing (static) tem- 
perature was found in tlie Sakai/Sugimoto model 33j. 
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lead to even solutions, which have the following expansion, 



1 



$4d(r)^ 8^ + 



(89) 



because only normalizable solutions are allowed, such that 
the constant leading term vanishes. The subleading coef- 
ficient M'^(r) has been multiplied by an additional factor 
— I/''$4d/8 for later convenience. Plugging above expansion 
into the eight-dimensional equation of motion, yields at lead- 
ing order in 



-9.(raO*4d = M'(r)$4d. 
r 



(90) 



This establishes that at least for a background geometry dual 
to a hydrodynamic expansion up to and including viscosity, 
the scalar meson equation is a Klein-Gordon equation with 
time-dependent mass. This might change when also encod- 
ing higher order effects like the relaxation time 19 , which 
currently appears to be out of reach of a supergravity ap- 
proximation 20 . 

We will now assess the error of the WKB approximation by 
plugging our meson solutions into the four dimensional Klein- 
Gordon equation. The error should be smaller than to be 
subleading to the viscosity contribution. We first determine 
the four dimensional meson solution by 



lim p^S<t){p,r)To[uj'-'^^] 



(91) 



where on the right hand side we plug in the mass ( 83 ) of the 
lowest holographic scalar meson solution, i.e., we set LU4,d ~ 

With (91) the error estimate Auj can be obtained from the 
Klein-Gordon equation 



1 



drTdr$4d{T) = {u;^'^\r) + Ac^ (t) ) " $4d (t) , (92) 



by linearizing in Acj. This yields 



£0 



-10/3 



+ . 



(93) 



which is sufficiently small, that is subleading to the viscosity 
terms arising from the geometry. (Also note that the fre- 
quencies Lj''^'(r) and the meson mass obtained from the holo- 
graphic expansion (89) agree up to and including order r~^.) 
However when encoding hydrodynamic effects in the geom- 
etry that are of sufficiently high order, we would be forced 
to consider a better approximation for our ansatz, e.g., by 
using higher order WKB. Moreover, beyond a certain order 
the WKB ansatz is expected not to work anymore because 
the coefficients S(f)i and S<f) are not expected to coincide to 
arbitrary order. 



IV. REGULARIZED DT ACTION 

In the static case, the free energy density can be related to 
the regularized D7 action by means of a Wick rotation. In the 
conventions of 13 the time direction is periodically identified 
with l3 = 1/T, such that the free energy density F is given by 




FIG. 7: The plot shows the relative difference of the masses 
of (type Uy) vector and scalar mesons for in — Sq 



with T the temperature and Vx ~ J d^x the (infinite) spatial 
volume of the boundary. The authors of 13 obtain 



NcNfT^X 
32 



12 



2mql\f\ 



' 6114m4 



y7 

(95) 



in the limit of low temperature. While it is clear that for 
a time-dependent background neither 1/T = J dt nor (|94| 
really make sense, we are still interested in how our action 
relates to above asymptotic result. 

When calculating the action Sd7, the integration along the 
holographic direction p on the brane is encumbered by UV 
divergences. The standard procedure of holographic renor- 
malization for the case of a probe D7 brane has been worked 
out in 31 . It consists in regularization by introducing a cut 
off pmax and addition of suitable counterterms formed from 
the induced metric and corresponding curvature on the slice 



P — Pmax . 



Pmax 

Sreg ~ J jO D7 dp + Li{pmax) 



(96) 



To be more precise, the procedure is usually formulated in 
terms of the coordinate 



^ = (p^+<I.^)-^/^ 

where the counterterms Li read 
1 



(97) 



Lb 
Lf 



48 



/7-R7 



ln(^™.„)*(a.7""^/79r + ^^/7^7)* 



(98) 



(94) 



I// is a finite counter term, with an arbitrary parameter a, 
that corresponds to different renormalization schemes. It can 
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be fixed in supersymmetric settings by the requirement tiiat 
the free energy vanishes. 7 is the induced metric on the z = 
Zmin siice and 



which agrees with ||95| 



'I' = arcsin 



(99) 



is the embedding coordinate of the D7 expressed as an angle 
on the internal 5*^. 

We express the counterterms in terms of pmax — 
$^(pmaa:))"^^^ usiug equatious (99) and (97) 

^ = -^rptnax + rm^^-i^^ + 0{p~^i^) (100) 

As a check we turn this back into z coordinates by self- 
consistently iterating (97) 



Pmax 



\ 



+ $2(- 



and using the D7 embedding (38). We obtain 

= ,2 , m^5+12«) 



4m c + 



12 



(101) 



(102) 



which up to the r-dependence arising from our rapidity coor- 
dinates is exactly the result of 13 , eq. (4.21), for a = —5/12. 
The renormalized action is thus 

Sren = N / dr Um / (Cdbi — P^t) dp 



+ ^m''(5 + 12Q)r, (103) 

where the last two terms are the divergent part of the coun- 
terterm, suitably rewritten as an integral, and the remaining 
finite part. Of course due to the renormalization, the limit 
can actually be carried out. For late times we obtain 



V. CONCLUSIONS 

The main goal of this article was to study fundamental 
fields in the holographic dual of an expanding viscous fluid. 
The dual geometry has strong similarity to a black hole geom- 
etry with time-dependent temperature. Small temperatures 
correspond to late-times. We determined the D7 embedding 
and calculated the consequences of dynamical temperature for 
the chiral condensate to three orders. The leading order gives 
the supersymmetric solution, the subleading corresponds to 
the adiabatic approximation and the subsubleading order in- 
cludes viscosity corrections going beyond the adiabatic ap- 
proximation. 

Moreover we calculated the meson spectrum and find that 
it agrees in the subleading order, though only the scalar 
mesons agree in the subsubleading order with the adiabatic 
approximation. The agreement crucially depends on the 
choice of ansatz defining the frequencies. We demonstrate 
that for our WKB ansatz solves the Klein-Gordon up to an 
error smaller than that inescapably introduced in the late- 
time expansion of the geometry. 

However, for the next order, which introduces the relax- 
ation time 19 , higher order terms would be needed in the 
WKB approximation in order to still keep that error smaller 
than the geometry's. 

It would be interesting to consider analogous properties for 
gauge theories which exhibit chiral symmetry breaking and 
hence are more closer to QCD. However up till now there 
is no description of an expanding plasma system in such a 
theory. 



Sren = N / dr T 



— ('12 I c^ I ^0 

(104) 

such that for a = —5/12 the configuration indeed relaxes to 
the supersymmetric setting. 

The above is our result. For comparison with (95) we per- 
form an additional, somewhat ill-defined step, namely replac- 
ing JdTT l/T. With L'^ = 2X£i and equation (34) we 
obtain for the supersymmetric scheme 



S-re 



3 • 21 



(105) 
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APPENDIX A: MESON SOLUTIONS 
(Pseudo-)scalar mesons 



5^0 



(8m*+9/[)^m^+3/[)*)£oi'® (l3m,*+12p^m^+3p*)eg' Tjoi* 1 



Wo - —For- + - 



5*1 = 



L2 

^o[t^i] 

2\/6m 



5*2 = ^o[W2] 



U)2 = 



4^/3r^ 



(-12m'' + 15p^m^+20p''m^+5/3°)eo-£'* (l9m^ -aS/J^m* -35p*m^ -Sp** )£o^*7)oi* p2_^2 
14m4(m2+p2)*^4/3 ' 14m* (m2 +p2)4^2 I" (m2+p2)2 

/_W3L!r2^\ ^ /i 2^0 \ 

(26m*-140p^m'^-22p''m*'+55p''m^+llp*)£oi'* 
30m*(m2+p2)5^4/3 

(41m*-269p^m^+121p^m,^ + llp*)eo^*7)oi^ 

30m'»(m2+p2)^T-2 

+ I, 6v/3m3 J \ r^'^^J ' 



4 o 2 2,4 
I m —op m +p 

(m2+p2)3 



5<I>3 = ^o[W3] 



(-134m^''+1506p^m*-1660p*m^-1045p^m''+342p*m,^+B7p^'')£oi'* 

154m4(m2+p2)«5T4/3 

(211m^°-2784p2m*+3585p*m'' + 1805p'^m''-912p^m^-57p^°)eo^''770-L* 
164m*(m2+p2)'iT-2 

I — m +op m —op m +p 
"I 7 — r, , o\4 



(m2+p2)'' 



4V^7T 



\ 77m3 J ^ ^2/3 . 



5*4 = J^o[oJ4,] 



(68ml^-1279p^ml°+3643p*m,*-630p''m8-1566p*m*+203pl''m^+29pl^)eoi* 
78m*(m2+p2)^r'*/3 

(l07ml^-2326p^m^''+7067p''m*-1840p^m^-3161p®m*+638pl''m^+29pl^)eo'''77oi* 

78m4(m2+p2)''T2 



m°-10p^m°+20p'^m*-10p'^m^+p 
(m2+p2)5 

29,/fi«r-V3eo' 



W4 = 



2V30n 
1,2 



13m,3 



1 2ao \ 

r2/3^; 



Vector meson (Type 112,3) 



r(22m''+21p2m2+7p*)£oi** (l3m'' + 12p2m2+3p*)£^/''77oi* , i 

(5a„2)o = ^o[a;o] I , „,^3,/,3 - 10^4(^2+^2)3,2 + 



Wo = 



2v^i7 
1,2 



30m4(m2+p2)3,4/3 
-4/3 



\ 16% 



v'2m3 



7t2/3 4/5^7 



(5aj/2)i = ^o[wi] 



OJl 



(-118m'' + 137p^m* + 164p*m^+41p'')eoi* 
126m*(m2+p2)*,r4/3 

(l9m^-3Bp^m''-35p''m^-5p^)eQ^*?70-£,* p2-m2 

14m4(m2+p2)''r2 (m2+p2)^ 



2V6T1 



41LfiT-*/3i 
21\/6m3 



41t2/3 ' 



13 



{5ail2)2 = ^o[<^2] 



(l78m*-S80p^m''-106p'*m*+315p''m2+63p^)eoi'^ 
180m<l(m2+p2)^T4/3 

_ (41m»-269p^m« + 121p«m2+llp»)e;^/*TOL** ml-Sp^m^+p* 
30m4(m2+p2)^T2 (m2+p2)^ 



(50/12)3 = •^o[W3] 



(-4666ml'' +48234p2m'*-52030p*'m''-29975p''m*'+9978p*m2 + 1663pl'')eo-£'^ 
4620m<l(m2+p2)6T'4/3 

(211ml°-2784p^m*+3585p*m'= + 1805p''m''-912p'*m^-57pl°)e|^^*'r7oi* 
154m4(m2+p2)''T2 



+ 



6,^2 4 ^4 2, 6 
-m +op m —op m +p 



(„2+p2)1 

4V5m , / 1663L°r-l/3eo \ /, 3420,,o \ 

- -L2- + 462y5m3 ^ ^ " 1663r2/3 4^ ^ 



(50/72)4 = ^o[<JJ4] 



(ll94ml^-20697p^ml''+55709p*m*~10890p''m''-22928p*m*'+2989pl''m^+427pl^)eo/'^ 

1170m-4(m2+p2)''r<l/3 

(l07ml^-2326p^mlf +7057p*'m*-1840p''m'5-3161p'*m''+638pl''m^+29pl^)so^*'>7oL* 

78m<l(m2+p2)^T2 

m^-10p^m'^+20p'*m^-10p°TTi^+p^ 



+ 



2V30m I 
1^4 = h 



(m2+p2)- 

/ 427L°r^'l/3eo \ ^ _ 8 70.70 ^ 
V 39x/30m3 y V 427r2/3 ) 



Vector meson (Type Ily) 

Vector mesons of type IIj, obey a different equation of motion (52), such that the WKB approximation yields a different 
result 



(Al) 



The corresponding holographic ansatz 

5 Ay — Saiiy T Ji{ J Lu dr) + Saiiy r Yi (J tj dr) 
gives rise to the following solutions for the five lowest mesons: 

(22m''+21p^m2+7p'')eoi* (llm'' + 4p2m^+p*)eQ''*J7oi 



2V2- 



L2 



+ 



30m<'(m2+p2)3^4/3 10m4(m2+p2)3 
^ - 7x2/™^) ' 



+ 



7L6T-^/3eo 
15V^m3 



{Saiiy)i = rJ"i[tJi] 



(-118m'' + 137p^m'' + 164p'*m^+41p'>)eo-t* 
126m-l(m2+p2)*T-4/3 

(81m<'-105p^m*-77p*m2-llp'')e„''*i70i'* p2_m^ 

42m4(m2+p2)''T-2 (m2+p2)2 



^ + 1 X (1 



1-2/3 4^y 



{Saiiy)2 = rJ"i[a;2] 



21^6.7 

(l78m**-880p2m''-106p*m*+315p'^m2+63p**)co/>*' 
180m'l(m2+p2)-V«l/3 

(l29m**-621p2m'=+40p''m''+209p''m2 + 19p'^)£o^*'.7oi* 
60m<l(m2+p2)^T2 



4 o 2 2,4 
I m — 3p m +p 

"1 ? — o . on3 



(m2+p2)3 



U!2 ~ 



4V3m. , 
—r^ 1" 



\ 10m3 J 



X 1 - 



38r,o 



21r2/3 4/5^7 ' 



(A2) 
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(-4666ml°+48234p^m^-52030p''m''-29975p'^m'*+9978p*m^ + 1663p^'')eo-t* 
4620m4(m2+p2)«r<l/3 

(3449ml°-34136p2m**+40675p''m'' + 15535p'>m*-8368p'*m2-523pl°)eQ'''770-L'' 
1540m<l(m2+p2)''T2 

_,_ -m''+6p2m''-6p^m2+p'3 



(m2+p2)* 

_ 4^m , / 1663L°T-^/3eo \ „ /i 31387)o \ 

^3 ^ V 462^m3 ) \ 1663r2/3 ' 



{5aiiy)i = rJ"i[a;4] 
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